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Abstract 

Signature change at high density has been obtained as a possible consequence 
of deformed space-time structures in models of loop quantum gravity. This article 
provides a conceptual discussion of implications for cosmological scenarios, based on 
an application of mathematical results for mixed-type partial differential equations 
(the Tricomi problem). While the effective equations from which signature change 
has been derived are shown to be locally regular and therefore reliable, the underlying 
theory of loop quantum gravity may face several global problems in its semiclassical 
solutions. 


1 Introduction 


Anomaly-free effective models of loop quantum gravity, derived for spherically symmetric 
configurations nua and cosmological perturbations at high density Pi, have revealed 
an unexpected phenomenon: At large curvature, signature change appears to be a generic 
feature of quantum space-time geometry as provided by this theory j5j. Not only the 
general phenomenon but also the specific form of signature change seems to be universal 
in these models, giving further support of the genericness of the effect. For the typical 
form of “holonomy modifications” used widely in models of loop quantum gravity, the 
speed of physical modes differs from the classical speed of light by a factor of (3(h,K ) : = 
a{h) cos(2AR) where a[h) > 0 is a function of the spatial geometry (the metric h), 6 
is a quantization parameter often assumed to be related to the Planck length, and K 
is a measure for extrinsic curvature (or the Hubble parameter in cosmological models). 
For large curvature K > 7t/(4<5), the speed is negative and commonly hyperbolic mode 
equations turn elliptic, which in these models are of the form 


d 2 (j) 

W 


P{h,K) 


A0 = S[(j)] 


( 1 ) 
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with source and lower-derivative terms S[(/>]. All modes </>, gravitational as well as matter, 
are affected in the same way. 

The overall picture shows some relationships with other approaches to quantum cosmol¬ 
ogy, mainly the no-boundary proposal of Hartle and Hawking [6] , and with other physical 
phenomena such as transonic flow or phases of nano-wires. Related mathematical ques¬ 
tions have been studied in the mathematical literature since the 1930s, following seminal 
work by Tricomi [Tj. Nevertheless, as a basis of effective space-time models, equations of 
the form (JT]) and the phenomenon of signature change show several new and surprising 
features. In this article, we give a general presentation of implications in cosmology. 

Section [2] gives a definition of signature change in the absence of a classical space-time 
metric, and contains a brief comparison with classical models of signature change. (A 
review of the formal origin and gauge independence of equation (jTJ) and signature change 
can be found in App. [A]) Section [3] provides several related examples in different areas 
of physics together with the appropriate mathematical formulation in terms of well-posed 
partial differential equations. Section 0 applies these results to gravitational questions, 
compares the pictures based on effective equations with wave-function methods, draws 
cosmological implications, and ends with cautious notes on global issues. 


2 Definitions 

Despite first appearance, the wave equation ([T]) is covariant under a deformed algebra 
replacing classical coordinate or Poincare transformations. The corresponding quantum 
space-time structure is not Riemannian, but has a well-defined canonical formulation using 
hypersurface deformations, as reviewed briefly in the appendix. 

A Riemannian manifold has a covariant metric tensor which transforms by g a 'b' — 
(dx a /dx a/ )(dx b /dx b ')g a b when coordinates are changed on the manifold. Accordingly, the 
line element ds 2 = g a bdx a dx b = g a iydx a dx b is invariant. Covariance under coordinate 
transformations (of solutions to the field equations) is canonically represented as gauge- 
invariance under hypersurface deformations in space-time [8]. Since the latter, along with 
Poincare transformations, are modified in the effective space-time structures we are con¬ 
sidering, they can no longer correspond to coordinate transformations. As a consequence, 
the effective “metric” does not give rise to an invariant line element. Instead of using 
metric-space notions such as geodesics and invariant scalar products of 4-vectors, in order 
to extract predictions there is only the possibility of computing canonical observables in¬ 
variant under the modified gauge transformations. In this way, one still has access to all 
observable information. However, the lack of a metric structure implies that many of the 
convenient and well-known techniques of evaluating solutions of general relativity are no 
longer available. 

Although there is no standard notion of geodesics and light rays or null lines, it re¬ 
mains straightforward to associate a causal structure to the modified space-time structure 
underlying (JT[), as long as (3 > 0. Instead of computing null lines for a metric, we just use 
characteristics of the wave equation (JT]) , for instance for gravitational waves to be specific. 
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(In the presence of inverse-triad corrections, scalar modes may propagate at speeds differ¬ 
ent from tensor modes ISIS!-) A characteristic is then a hypersurface which at any point is 
normal to a vector k a satisfying [k 1 ) 2 — (3\k\ 2 = 0. With these characteristics instead of null 
lines, we can define light cones, a causal structure, and derived notions such as Penrose 
diagrams. 

Characteristics exist as long as (3 > 0. For holonomy modifications, /3 changes sign at 
high density, and the characteristic equation no longer has non-trivial solutions k a ^ 0. 
As a consequence, (CQ) does not provide a causal structure in such a regime. Since the 
same condition of /3 < 0 makes the mode equation (pQ) turn elliptic, we call this regime 
“Euclidean,” while for /3 > 0 we call it “Lorentzian.” In this way, we generalize the two 
standard discrete choices of signature to a continuous range of values taken by /? as it 
varies from the classical limit (3 = 1 at low curvature to a negative value at high curvature. 
Unless (3 = ±1, the canonical fields from which (JT|) is derived do not provide a classical 
notion of Lorentzian space-time or 4-dimensional Euclidean space. But the most important 
properties regarding physical consequences, including the existence of a causal structure 
and the type of initial or boundary value problem required for reliable solutions, only 
depend on the sign of /3 rather than its precise value. For this reason, we still speak 
of Euclidean or Lorentzian signature even if we have neither Euclidean nor Lorentzian 
space(-time). 

Signature change has been studied in quite some detail in classical general relativity 
(/l = 1); see for instance [XU] - However, such models are crucially different from what is 
considered here, in that they have /3sgn(det h ab ) (or /3sgmA 2 ) changing discontinuously. As 
a consequence, in classical models the transition from Euclidean to Lorentzian signature is 
always singular, which is not the case in our effective models. (Note that for this reason 
signature change in the model of na is not analogous to deformed space-time structures.) 
The subtleties and controversies related to the distributional nature of solutions with clas¬ 
sical signature change, discussed for instance in [T2| [13.; XU QA, T6[ [TTj '18], do not play a 
role in our context. (The models we study could be considered as a version of dynamical 
signature change as anticipated in the conclusions of [T9].) 

It is necessary to consider inhomogeneity in order to see modified space-time struc¬ 
tures in which modes obey (pQ). In particular the phenomenon of signature change, which 
is the main topic of this article and is realized because /3 may change sign, can only be 
seen in inhomogeneous models: it is manifested by the relative sign in space and time 
derivatives in field equations. However, signature change (or a modified space-time struc¬ 
ture) is not a consequence of inhomogeneity, which in the perturbative context would be 
dubious [2D]]. The origin of signature change lies in the modification that one makes in 
the classical dynamics if one quantizes the theory following loop quantum gravity, which 
provides operators for holonomies (exponentiated and integrated connections) instead of 
ordinary connection components. This modification appears already in the background 
dynamics of a homogeneous minisuperspace model, but in this context one does not notice 
signature change because all spatial derivatives vanish. Perturbative inhomogeneity then 
makes the effect visible. Still, inhomogeneity or perturbation theory is not the origin of 
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signature change, which is easy to see by the following two arguments: First, signature 
change persists in a perturbative field equation no matter how small the inhomogeneity is, 
as long as it is non-zero. The coefficients of space and time derivatives in (JTJ) depend on 
background quantities, and for ft < 0 they have the same sign for all values of inhomogene¬ 
ity. Secondly, the same form of signature change appears in spherically symmetric models 
in which inhomogeneity need not be treated perturbatively. 

3 Related effects 

In the explicit form as provided by models of loop quantum gravity, signature change is a 
new effect. But it has precursors in physics as well as mathematics. Besides the examples 
discussed in details below, the case of helically symmetric binary systems [21] is worth 
mentioning. The configuration is described by partial differential equations which change 
signature along a spacelikc direction (far from the center), rather than a timelike one as in 
our cosmological models. 

3.1 Hart le—Hawking wave function 

A Lorentzian space-time cannot be closed off at a finite time without a boundary. As 
proposed in [6J, however, one may postulate that quantum gravity gives rise to a mod¬ 
ified space-time structure that allows a transition to Euclidean 4-dimensional space. A 
Euclidean cap can then be attached to a space-time manifold which becomes Lorentzian 
at low curvature. In [6] and the literature based on it, this scenario has been used mainly 
to specify the high-density asymptotics of wave functions satisfying the Wheeler-DeWitt 
equation of homogeneous models. 

The scenario suggested by effective constraints in models of loop quantum gravity can be 
seen as a concrete realization of the quantum-gravity effects that may give rise to signature 
change. The condition that extrinsic curvature vanish at the interface between Euclidean 
and Lorentzian parts of semiclassical solutions 122 is then replaced by /3(h, K) = 0. Never¬ 
theless, it is not guaranteed that the same consequences are implied for wave functions, not 
the least because the possibility of signature change relies on a quantum-geometry effect 
(based on the use of holonomies in loop quantum gravity) which turns the Wheeler-DeWitt 
equation into a difference equation [23, [24]. The minisuperspace discreteness implied by 
this difference equation is relevant especially at high density, that is in the Euclidean phase 
where asymptotic properties are discussed according to [5]. Qualitatively different conclu¬ 
sions for wave functions could then be reached, so that it is not clear whether the close 
relationship in the space-time picture of [6] with the models discussed here implies a similar 
relationship in predictions. Specific details of wave functions might well be closer to [25] 
than to [6]. We leave this question open in the present article, in which we are concerned 
mainly with effective equations. 
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3.2 Analog condensed-matter models 


In addition to related cosmological models in terms of space-time properties, there are 
rather different physical phenomena which give rise to mathematical descriptions similar to 
(JIT). The main examples known to us are transonic flow and hyperbolic metamaterials [26]. 
Especially the latter phenomenon provides an interesting analog picture of the cosmological 
effects. 

An example of a hyperbolic metamaterial is given by an array of parallel conducting 
nano-wires immersed in a dielectric medium. Such a material behaves as a conductor 
in the direction parallel to the nano-wires while in the normal one it exhibits dielectric 
properties [27]. Choosing the £-axis to be parallel to the nano-wires, the wave equation for 
the component E z =: 0 of the electromagnetic held can be written as [27J 

d 2 (j) 1 d 2 (j) 1 / d z (j) d 2 0\ = 

dt 2 e± dz 2 (|| \dx 2 dy 2 ) 


For the “wired” metamaterials the electric permittivity is e\\ > 0, as for the standard 
dielectric medium. However, because the metamaterial is conducting in the ^-direction, e± 
will be negative in a sufficiently low frequency range!)] In particular, for the Drude theory 
of a conducting medium 


e± = 1- 


CJp 

UT 2 


(3) 


where cup is the plasma frequency dependent on the geometry and composition of the 
metamaterial. At low frequencies (c o < cup), e_i_ > 0 and the spatial part of equation ([2]) is 
elliptic. However, for oj > cap, we have e± < 0 and the sign in front of the second derivative 
with respect to z changes. The spatial Laplace operator becomes hyperbolic, with the 
^-component playing a role of the second time variable. The dispersion relation 



e± £|| 


0, 


(4) 


where k jj := k 2 -\-k 2 , for fixed co is in this regime no longer represented by ellipsoids but 
hyperboloids. In consequence, the k -space undergoes a topology change. 

Passing between the elliptic and hyperbolic regions is not only a matter of frequency 
dependence. For a fixed (sufficiently low) frequency, the transition between the different 
signs of 6j_ can be induced by a rearrangement in the structural form of the “wired” meta¬ 
material. In particular, melting of the nano-wires to the liquid phase (associated with a 
first-order phase transition) has been observed experimentally. During the process the value 
of ej_ smoothly changes its sign (from negative to positive) with increasing temperature. 

assume here that the magnetic properties of the metamaterial are as usual, that is magnetic 
permeabilities in all directions are positive. However, it is worth mentioning that metamaterials with both 
electric permittivity e < 0 and magnetic permeability y < 0 in some frequency band have been constructed. 
An interesting property of such metamaterials is that the refractive index is negative n = — y/ey, leading 
to very interesting and sometimes counterintuitive behavior [55]. The relevance of this phenomenon in the 
context of analog models of quantum gravity is an open issue. 
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Alternatively one can imagine the transition to be of the second order. This idea in the 
context of signature change due to holonomy corrections has been discussed in [29l 1301 l3Tj . 
Let us imagine that the nanowires are immersed in a dielectric fluid, having an ability to 
rotate. Furthermore, let us introduce magnetic dipole-dipole type interactions between the 
wires by attributing magnetic moments to them. Then, in the high temperature unordered 
phase, with all nano-wires pointing in random directions, there is hardly any electric con¬ 
ductivity. But when the phase changes to an ordered one (lowering the temperature), a 
significant electric current starts flowing, just as time starts flowing in our universe models 
when the density is small enough to trigger signature change to Lorentzian space-time. 

In the gravitational case, the process can be modelled by a bi-metric gravity theory, 
where the effective metric experienced by the fields is 

9nu $fii/ ■ (5) 

The metric may be viewed as a 4-dimensional analog of the inverse of the electric permit¬ 
tivity tensor e t] . The x^ is an effective mean field having the interpretation of an order 
parameter, defined such that in the unordered phase |x| = 0 (here |x| = \J <3^ XvXv) ■ 
In this case the fields experience the Euclidean metric g^ = 8^. In the fully ordered 
phase |x| = 1, the Lorentzian signature emerges in a spontaneously chosen direction in the 
4-dimensional space. 

In order to describe the transition from the unordered to the ordered state more quan¬ 
titatively let us postulate the following form of the free energy for a model with a test field 
0: 


F 


dV 5^ + 


- - 2|xP 


dn(pd v (t) + C 


— 



-i) IVI 2 + |lxl 4 

V 


where C is a constant. Because of the x^X v 'd^d^ contribution to the kinetic term, 
F is not SO(4)-invariant in the internal indices of yA The explicit symmetry-breaking 
term can be treated perturbatively. Then, equilibrium corresponds to a minimum of the 
SO(4)-invariant potential V(x^iP)- In order to model signature change in loop quantum 
cosmology with maximum energy density Pqg, the parameters of the potential are fixed 
such that at energy densities p > /?qg the vacuum state maintains the SO(4) symmetry 
- the order parameter is |x| = 0. For p < pqg, we choose the potential so that its 
minimum is located at |y| = ^/l — p/pqg, with some spontaneously chosen direction in 
the 4-dimensional configuration space of the order parameter. Without loss of generality, 
we may assume x to point in the 0-direction (xo = |x| and x% — 0 for i — 1, 2, 3); thus, 
goo — 1 — 2xoXo = —1 + 2p/pQQ = —(3 and ga = 1. The corresponding wave equation for 
the test field is 

-^ +A0 = 0 ’ (6) 
which is of the form ([]]) with (3 = 1 — 2p /pqq as in models of loop quantum cosmology. At 
|x| = 0 ((3 = 1) the equation exhibits the R A x SO(4) symmetry while in the fully ordered 
state |x| = 1 (/3 — —1) the A 1,3 x SO(l,3) Poincare symmetry is satisfied. 
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It is worth stressing that the presented model relating signature change with the spon¬ 
taneous symmetry braking of (SO (4) has no support from more fundamental considerations 
at the moment. Therefore, its physical validity requires further studies. In particular, a 
possible relation of the order parameter with the elementary degrees of freedom of loop 
quantum gravity remains unclear. Nevertheless, the model presents an interesting example 
because of its formal resemblance to signature change in loop quantum cosmology. There 
is, however, an important conceptual difference: while analog models of signature change 
are fully well-defined and even have some observable properties, quantum-cosmology ef¬ 
fects on space-time itself (rather than matter fields in space-time) may be subject to several 
global problems on which we will comment in Sec. 14.2.61 


3.3 Mixed-type characteristic problem and Tricomi equation 

Partial differential equations of mixed type have been studied in the mathematics literature 
since the 1930s, beginning with the Tricomi problem for the differential equation 


d 2 u d 2 u 
dy 2 ^ dx 2 


(7) 


in the (x, y)- plane. This type of equation can be seen as an approximation to our equation 
dU) around one of the finite boundaries of the elliptic regime, where (3 = 0. Here, we recall 
relevant features of a suitable characteristic problem mm- 

Let us suppose that we would like to solve equation (J7|) in a domain 12 = lb Ul2 2 , where 
12i and 122 are elliptic and hyperbolic domains, respectively. For the Tricomi problem, the 
sub-domain lb cannot be chosen in an arbitrary manner. Namely, the boundary <91b 
is determined by characteristics of (J7J) attached to the endpoints of <912i at the interface 
between the hyperbolic and elliptic regions; see Fig. [T] The characteristics in the hyperbolic 
region are determined by the characteristic equation to (JTJ): (dx) 2 + y(dy) 2 = 0. This 
equation has two families of solutions with a different sign, enabling us to introduce the 
characteristic coordinates 


« := z + ( 8 ) 

ri '■= x - |(-j/) 3/2 . (9) 

(We have y < 0 in the hyperbolic region.) 

Solutions to the Tricomi equation (J7J) are unique and stable if adequate boundary 
conditions are imposed. Applying the notation from Fig. [U they are as follows: in region 
lb, the value of u has to be specified at the boundary curve A, while in region 12 2 the value 
of u at the characteristic curve C has to be given. The value of u\b is a result of fixing the 
boundary condition u\cua- 

The reason why the boundary condition in the hyperbolic region is imposed on the 
characteristic curve becomes clear by analyzing the case of the standard wave equation 
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Figure 1: Graphical representation of the characteristic problem from the mixed-type 
Tricomi equation. Solutions to a partial differential equation (JTj) (hyperbolic for y < 0 
and elliptic for y > 0) exist uniquely and stably in the shaded region U ff 2 , provided 
boundary data are imposed on one characteristic C of the hyperbolic region together with 
a curve A in the elliptic region. Characteristic coordinates (HI]) are indicated by £ and 77 . 





with one spatial dimension: 


d 2 u d 2 u _ 
dy 2 dx 2 

The corresponding characteristic equation (dx ) 2 — (dy) 2 = 0 has solutions x±y 
allow us to introduce the characteristic (light cone) variables 

f : = x + y, 
rj := x — y 

In terms of these variables the wave equation (fTOl) reduces to the canonical form 

d 2 u 
d^dy 

with general solution u = /(£) + g(y). An important observation is that an initial-value 
problem (Cauchy problem) in the variables (x, y) translates into the boundary-value prob¬ 
lem in the characteristic variables (£, 77 ). In particular, the Cauchy-problem initial condi¬ 
tion u\ y= o and du/dy | 0 translates into the Dirichlet boundary condition at the light cone 

m|j = o and u\ v= q. 

Using characteristics, it is therefore easier to combine hyperbolic and elliptic regimes. 
In the case of the Tricomi equation, the boundary conditions in the elliptic domain have 
to be smoothly extended into the hyperbolic sector. This requires the introduction of 
characteristic variables in the hyperbolic sector, imposing boundary conditions at either 
£ = constant or y = constant (but not both). 

4 Cosmology 

Signature change plays an important role in cosmological scenarios. In practical terms, one 
of the most important consequences is the presence of instabilities, as indicated formally 
by an imaginary speed of sound in mode equations. The same kind of instability plays a 
role in the mathematical discussion of equations of the form (JT]) : One of the conditions for 
a well-posed initial/boundary-value problem of a partial differential equation is the stable 
dependence of solutions on the chosen data, in addition to the conditions that a solution 
exist and be unique for given data. 

4.1 Instability 

Stability is easily violated when one attempts to use an initial-value problem for an elliptic 
equation. For instance, Fourier modes of solutions to the standard wave equation would 
not oscillate as per exp(±icot) but change exponentially according to exp(±cut). The ex¬ 
ponentially growing mode implies instability in the sense that solutions, if they exist and 
are perhaps unique, depend sensitively on the initial values. Choosing a boundary-value 
problem in t (as well as spatial directions), on the other hand, ties down the solution at 


(13) 


( 11 ) 

( 12 ) 


( 10 ) 
= C, which 
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both ends of the f-range, so that the growing branch is sufficiently restricted for solutions 
to be stable. 

Signature change implies instability of initial-value problems, but it is stronger in two 
respects. First, signature change affects all modes (of matter or gravity) equally, which 
all become unstable at the same “time.” It is therefore a space-time effect, rather than 
an exotic matter phenomenon. Secondly, it can be seen in space-time symmetries such as 
hypersurface deformations or the Poincare algebra. None of these effects happen in known 
cases of instabilities in matter theories or higher-curvature gravitational actions, such as 
[33] . Signature change is therefore more fundamental than other phenomena hat might 
imply instability. An important property is the fact that the theory does not provide any 
causal structure whatsoever when the signature turns Euclidean. 

4.2 Scenarios 

When used in cosmological model-building, signature change implies a new and interesting 
mixture of linear and cyclic models. One of its main consequences can be described as a 
finite beginning of the universe. In the simplest versions of existing models of loop quantum 
cosmology, effective equations do not imply divergences at high density but instead extend 
solutions to a new low-density regime [34]. The main example is a modified Friedmann 
equation 

sin(h'H ) 2 87 tG 

I 2 = ~1T P 

postulated for spatially flat isotropic models sourced by a free, massless scalar 0. The 
energy density p is therefore of the form p = -(p^/a 6 with the momentum p^ of 0 and the 
scale factor a. In (fT4]h 6 is a parameter which approaches 5 —> 0 in the classical limit, 
and whose precise value remains undetermined. (As a parameter motivated by quantum- 
gravity effects, it is often assumed to be close to the Planck length.) Moreover, "H = 
(2<5 )“ 1 arcsin(2<5a/a) is a modified version of the Hubble parameter. It is straightforward 
to see that the modified equation (fT4|) implies that the energy density, for fixed 5, is always 
bounded, p < Pqq := 3 /( 87 tG5 2 ). The parameter pqq is the maximum kinetic energy 
density realized at the bounce. There is then a good chance that cosmological singularities 
may be resolved. Indeed, the quantum model with a free, massless scalar can be solved 
completely [35], implying that the scale factor 

a(t) = a 0 (l + 247 tGpq G t 2 ) 1/6 

never becomes zero. (See also [36].) 

Equation (fT4|h although it remains unclear how generically it is realized for the back¬ 
ground evolution [37] of general homogeneous models, has suggested a picture in which 
our expanding universe descends from a preceding collapse phase which produced large 
but not infinite density at the big bang. However, when modes on such a bouncing back¬ 
ground are subject to equations of the form ([I]), as required for a consistent system of both 
background variables (a, 0 ) and inhomogeneous modes in models of loop quantum gravity, 


(15) 
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the transition is not deterministic owing to the lack of causal structure at high density. In 
the cosmological case, the role of K in ([I]) is played by the modified Hubble parameter 7~L. 
With P(l-L) oc cos(25"H), one can easily see that signature change happens when the energy 
density is half the maximum it can achieve (or half the bounce density of a homogeneous 
model). Given the solution (fT5lh the mode equation ([I]) is elliptic for 


247tGpqg 


(16) 


(For pqg close to the Planck density, the Grange in which the equation is elliptic is close 
to the Planck time, but it can be significantly longer in models with a scalar potential 
or strong quantum fluctuations at high density [38].) The high-density regime cannot 
be accessed causally, and there is no deterministic transition from collapse to expansion. 
For practical purposes, such a scenario therefore shares some features with the traditional 
singular big-bang model, in which one must pose initial values just after the singularity. 
Similarly, one must pose initial values “after” the Euclidean phase of a signature-change 
scenario. However, the relationship between singular and signature-change models is subtle, 
as shown by the detailed discussion of suitable data for mixed-type differential equations 
on which we now embark. 


4.2.1 Gravitational waves 

To be more specific, let us focus now on the case of gravitational waves with holonomy 
corrections. In this case, the equation of motion for each component of polarization of the 
gravitational waves is [39j 


d 2 </> 

3t 2 



90 P A , n 

— -—A0 = 0 , 
at a 2 


(17) 


where H := ^ is the Hubble parameter and the deformation factor is 

P = cos(2 6U) = 1 - 2— . (18) 

Pqg 

Depending on the value of the parameter /?, equation ora can be either hyperbolic (/? > 0), 
elliptic (P < 0) or parabolic (P = 0). The transition between elliptic and hyperbolic type 
(associated with signature change) takes place at p — \pqg- At this moment, the square 
of the Hubble factor H reaches its maximal value H 2 vax = ^ttGpqq. For the model with 
a free scalar field introduced in the previous section, this takes place at t — ±t max = 
±(247rGp QG )- 1/2 . 

Equation (1X71) is precisely of the form ([1]) with the source term 

S M = -(3J7-f)* (19) 


11 





Because generally $ ^ 0 at (3 = 0, signature change is associated with a divergence of the 
source. The divergence, however, does not lead to pathological behavior at the level of 
solutions to equation because, as we will see later, it is due to a regular pole, which 
does not disturb regularity of the solution. In the vicinity of (3 = 0, equation (fTTll can be 
approximated by 


d 2 0 

dt 2 




( 20 ) 


which after double integration leads to the solution 


(j) = Ci + c 2 



( 21 ) 


where C\ and C 2 are some functions of the spatial variables. Because f3 occurs in the 
numerator of the integrand, the approximate solution is indeed regular. 

Let us now analyze the behavior of equation (1X71) in the vicinity of the instance of 
signature change for the model with a free scalar field. Without loss of generality, we 
perform an expansion around t = — t max ? for which 



- 5 tjj 
t 4 -1 4 


t + Li 


+ 


t +t n 


4 t 2 


+ O ((t + tmax) 2 ) • (22) 


(In the first step we have used (fl5lh (1T61) and (fTSl) in order to write a(t ) = ao(l+ t 2 /tm ax ) 1//6 
and /3(t) = 1 - 2(1 + 1 2 /1^)- 1 .) Because (t, + f max ) (sH - f3/f3^J = -1 + O (t + t max ), the 
pole at t = —tmax is of the regular type. Furthermore 




(t/tmax) 2 - 1 

(1 + (t/tmax) 2 ) 4 / 3 


t + t n 


2V3f 

^ t/ri 


+ O ((t + tmax) 2 ) , 


(23) 


where we fixed a 0 — 1. 

Applying the above expansions to the leading order and reducing to the (l+l)-dimensional 
castlE we obtain 

d 2 </> 1 d<j) t + t max d 2 (p _ 

dt 2 t + t max dt 2 1 / 3 t max dx 2 

Redefining the time variable to 

!,:= (2V3 i L) 1/3(( + (m " ) ’ <28) 

equation (1241) can be rewritten as 


d 2 (j) d 2 4> 1 dcj) 

dy 2 ^ dx 2 y dy ' 


(26) 


2 Of course this reduction is formal only because gravitational waves do not occur in (l+l)-dimensions. 

Alternatively, one can consider the reduction as being a result of introducing translational symmetry in 

the y and z directions, leading to 4>(t, x, y , z ) = </>(f, x). 
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The left-hand side of this equation is identical to Tricomi’s expression in (J7J) while the right- 
hand side can be treated as a source term. Applying the characteristic Tricomi problem 
to obtain solutions to the equation (l26l) is therefore justified. 

It is worth noting at this point that equation (12611 in the vicinity of y = 0 reduces to 


d 2 (j) 1 dcj) 

dy 2 ydy 


(27) 


with solution </> = C 1 +C 2?/ 2 (in agreement with (|2ip ). The solution is perfectly regular across 
the signature change, in contrast to classical signature change. In particular, signature 
change resulting from the line element ds 2 = — tdt 2 + dx a dx“ has been widely studied in 
the literature ixa using- In this case, the analog of equation (127)) is = w ^h 

the metric signature change at y — 0. While both forms of the equation look very similar, 
the extra factor of “2” plays a significant role leading to the non-differentiability of solution 
at the signature change. For discussion of controversies around this issue as well as some 
proposals of dealing with the problem we refer to Ref. [I6lj . 

Coming back to expression (1261) . the corresponding equation for the Fourier component 
4>k = (27 r)~ 1//2 f d ke~ lkx 4>{x) is 


d 2 4>k 
dz 2 


1 dcpk 

z dz 


- Z(j) k 


0, 


(28) 


where z := y\k : J> has been introduced to absorb the wave number k. In order to facilitate 
solving this equation, we observe that by differentiating the Airy equation 


d 2 u 
dz 2 


— zu = 0 , 


and subsequently substituting for u using again (129|) . we find an equation 


(29) 


d 3 u 1 d 2 w d u 

d z 3 z dz 2 dz 


(30) 


identical to (!28l) if (j) k = du/dz for a suitable (fc-dependent) u solving (|29|) . Solutions to 
(I28p are therefore 

4> k = A k Pa(z) + B k Bi'(z ), (31) 

where Ai'(z) and Bi'^) are the Airy prime functions (derivatives of the Airy functions) 
and A k and B k are fc-dependent constants of integration. Because of the reality condition 
for the field 0 the following relations have to be fulfilled: A* k = A_ k and Bl = B_ k . 

As an example, we will now apply the Tricomi boundary conditions to the resulting 
2 -dimensional solution 


/ + °° Ah 

-j=e kr - [A t M (y\kn + B k Bi' (j/lil 2 ' 3 )] . (32) 

-00 V 
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First, let us choose the contour A : = {(x, y) G R 2 , y = 0, — tt < x < i r} and fix 0|^ = cosx. 
This condition reduces (l32|) to 


c f>(y , x ) = cos(x) 


1 — c 0 Bi'(0) 

——Ai (j/) + c 0 Bi (j/) 


(33) 


with a parameter cq G R. Secondly, the form of A implies that: 


C := <( (x,y) G R 2 ,x = vr - -(-y) 3/2 , 


3tt\ 

~2 J 


2/3 


<y <o >, 


b ■= 


Finally, by choosing 


(x,y) G M 2 ,a: = -tt + jj(-//) 3/2 ,- < V < 0 


01 C = COS ^7T - ^-y) 3/2 ^j Ai '(2/)/ Ai, (°) 
the constant Co is hxed to equal zero, which leads to 

(j>(y,x) = — 3 1 ^ 3 r(l/3) cos(x)Ai , (//). 


(34) 

(35) 

(36) 

(37) 


The solution together with the boundary d£l = A U B U C is shown in Fig. [2j For all 
other values of c 0 , the contribution from Bi \y) in (l33l) does not vanish, and therefore the 
solution is characterized by an exponential growth in the y > 0 region. Nevertheless, in 
the domain enclosed by dVt = AU B U C the solution remains regular. It is transparent 
from Fig. [2] that the monotonic solution in the elliptic region (y > 0) transforms smoothly 
into the oscillatory solution in the hyperbolic domain (y < 0), as expected intuitively. 


4.2.2 Characteristic problem 

We are now in a position to adapt the Tricomi problem to effective space-time models of 
loop quantum gravity. We begin with a qualitative discussion which also provides heuristic 
reasons for why the mixed-type initial-boundary value problem must be used. 

In a non-singular universe model in which collapse is joined by expansion, one can start 
with initial values in the collapsing Lorentzian phase in the past and evolve all the way 
to the boundary with the Euclidean phase (at t = — f max in the solvable model). For a 
second-order hyperbolic differential equation in this regime, we need to fix initial values 0 
and (j)', where we denote by 0' the derivative of 0 normal to constant-f hypersurfaces. In 
this regime, 0' may therefore be considered as the time derivative of 0. When we reach 
the Euclidean phase, time no longer exists in a causal sense and our coordinate t is one of 
the four spatial ones. We will continue denoting the f-derivative as 0', now strictly in the 
sense of a normal derivative replacing the time derivative. 

Once we reach the Euclidean phase and enter the range between t = — f max and t = f max 
in the model, we must switch to a boundary-value problem for wcll-posedness. We can 
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Figure 2: Density plot of the solution (1371) together with the boundary dfl = AUBUC. This 
special solution is exponentially decreasing for y > 0, as can be seen by the fading-away of 
oscillations. 
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Figure 3: Illustration of indeterministic behavior implied by the mixed-type differential 
equations of signature-change cosmology. The boxed data for a held 0 — the held and 
its hrst normal derivative 0' in the past Lorentzian phase as well as the held 0 (but not 
its normal derivative) at the beginning of the expanding phase — must be prescribed. 
The other data then follow from a solution to the mixed-type differential equation: 0 and 
0' at the left boundary of the Euclidean region by an initial-value problem in the past 
Lorentzian phase (left arrow); 0' at the right boundary of the Euclidean phase as the 
limiting normal derivative of the solution in the Euclidean phase based on the evolved 0 
on the left boundary and the prescribed boundary data 0 on the right boundary (center 
circle); and finally the future held 0 and 0' in the expanding phase evolved from initial 
values at the right boundary of the Euclidean phase (right arrow). This initial/boundary- 
value formulation is made more precise in the Tricomi problem, shown in Fig. [H 
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use our field 0 evolved up to the part of the Euclidean boundary bordering the collapse 
phase, at f = t max . But we must complete it to a set of boundary values on a hypersurface 
enclosing the Euclidean region in which we are looking for a solution. For a solution 
in all of space(-time), this complete boundary includes asymptotic infinity in directions 
other than t as well as the boundary at t = f max where the Euclidean region borders on the 
expanding phase. We may take for granted that asymptotic fall-off conditions for our fields 
are imposed at spatial infinity, but we must choose a new and arbitrary function on the 
border with the expanding phase. When this function is fixed, together with the evolved 
past initial data, we obtain a solution in the Euclidean phase. We may then compute 
normal derivatives on all constant-f hypersurfaces in this region and take the limit toward 
the border with the expanding phase. We obtain suitable data for an initial-value problem 
in the latter. In summary, as illustrated in Fig. [3l in addition to past initial data we must 
specify one additional function for every mode at the beginning of the expansion phase. 
Past initial data in the collapse phase do not determine a unique solution in the expansion 
phase. There is no deterministic evolution across the Euclidean high-density regime. 

The intuitive description just given provides the correct picture of free and determined 
data, but it is not the most suitable one for the problem. One question one may pose is 
about the smoothness of solutions: A solution to (Op) should be required to be smooth (or at 
least twice differentiable) at the borders between Lorentzian and Euclidean regions. In the 
picture of Fig. |3] it is not clear whether this requirement is satisfied. We have to match two 
solutions, one in the past Lorentzian phase and the one in the Euclidean phase, so that the 
normal derivative 0' obtained from these solutions changes differentiably across the border. 
Since the field 0 on the border with the expanding phase is free and affects the solution 
in the Euclidean region but not in the collapse phase, differentiability seems in danger 
unless one severely restricts the boundary data left free around the Euclidean region: If 
we keep past initial data fixed but vary the boundary field 0 at the border t = t max with 
the expansion phase, the limit of 0' at t = —t max will change for solutions of the elliptic 
equation. Only a narrow range of boundary fields 0 would seem to provide a 0' matching 
the evolved initial data in the collapse phase. If this is so, one would essentially be led back 
to an initial-value problem across the whole domain, because 0 and 0' everywhere would 
be uniquely determined by initial values, a conclusion which cannot be correct if there is 
an elliptic regime. By these considerations, for a smooth solution we would have to give 
up some of the initial values in the collapse phase as free data, but using just 0 and letting 
0' be determined by smoothness cannot be correct either, because we would end up with 
a boundary-value problem which is not well-posed in the Lorentzian phase. 

In order to address differentiability, a Tricomi-style characteristic problem is more ap¬ 
propriate. In the Lorentzian phase, one specifies data not on initial-value surfaces but on 
characteristics of the differential equation as long as it is hyperbolic. Far away from the 
border with the Euclidean region, these characteristics are standard light rays, but they 
approach the border normal to constant-f surfaces and end there. One example is shown as 
curve C in Fig. SI Considering just one spatial dimension, a second characteristic starting 
from the same point as C completes the first one to a characteristic triangle with the two 
characteristics as well as the enclosed border with the Euclidean region as its sides. 
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Figure 4: Well-posed characteristic problem for a mixed-type partial differential equation: 
The function is specified on one characteristic arc C of the equation in the Lorentzian 
phase and on an arc A in the Euclidean phase connecting the endpoints of C and of 
a second characteristic arc starting at the same point as C (dashed). The data on the 
second characteristic arc are determined for solutions of (JT]) and cannot be chosen freely. 


The problem is well-posed [?, [32] if one specifies the field </> on one of the two char¬ 
acteristics and on an arc A which starts and ends at the border-endpoints of the two 
characteristics but stays in the Euclidean region in between. If these data are smooth 
on the curve obtained as the union of C and A one obtains a unique, stable and smooth 
solution in the interior of the union of the characteristic triangle and the region enclosed 
by the curve A. Smoothness is guaranteed in this way, with data on the curve A being 
completely free except that they must extend the data on C smoothly. The indeterministic 
behavior (requiring future data) is therefore confirmed. If one lets A approach the other 
border of the Euclidean region, one can proceed as in the intuitive description based on 
Fig. [3] An initial-value problem starting at the beginning of the expanding phase is then 
well-posed toward the future. 

In four dimensions, the characteristic triangle is extended to the interior of the future 
light cone starting somewhere in the collapse phase, and the region enclosed by a curve A 
is extended to a dome enclosed by a 3-dimensional hypersurface in the Euclidean region. 
Although we are not aware of a complete mathematical proof, one may expect the charac¬ 
teristic problem to be well-posed with data given by the mode on half the light cone and 
on the hypersurface in the Euclidean region. 

The smoothness of solutions shows that there is still a well-defined manifold which 
combines the two Lorentzian phases and the Euclidean phase, but this manifold is no 
longer Riemannian. Standard obstacles to signature change in general relativity therefore 
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do not apply. The use of a finite region in the characteristic formulation does not cause 
a problem in cosmology because we always have observational access to a finite part of 
our universe. The characteristic triangle can always be chosen large enough to include the 
whole observational range. 

4.2.3 Boundary conditions at the hypersurface of signature change 

Besides the boundary, initial-value and the characteristic problems discussed in the pre¬ 
vious subsection we would like to present one more possibility of posing the initial-value 
problem for the evolution in the Lorentzian domain. Let us start by finding the solution 
in the Euclidean domain. 

In this domain (which we call Qe) the equations of motion are of the elliptic type. 
Therefore, for finding solutions in a subset C VLe the proper boundary conditions have 
to be imposed. Specifying 4>\a o (where dQ is the boundary of f2) is sufficient to determine 
a solution <f> E to the elliptic equation in the domain uniquely. An interesting situation 
corresponds to the case when the boundary dQ, is expanded so that it encloses the whole 
Euclidean domain: dQ, —>■ <9f Ie- Such a boundary can be decomposed as follows ()Qe = 
<9f2 + + + <9f2oo- The dtt± are boundaries between the Euclidean and Lorentzian 

domains at ±t max , respectively, while <9f2oo encloses the Euclidean domain at spatial infinity 
(assuming that the space part is unbounded). The solution (f>E at Qe obtained by imposing 
boundary (Dirichlet or von Neumann) conditions at 8£Ie can now be used to determine 
initial (or final) conditions for the Lorentzian sectors fl±, where “ + ” corresponds to 
t > +t max and “ — ” to t < — f max . The ±f max (as previously) are the times at which 
signature change is taking place. 

Let us decide to impose the Dirichlet boundary condition leading to a solution 

c f>E■ Then, the Cauchy initial conditions at ±f max for the solutions c/)l± br the Lorentzian 
domains Q l± respectively are: 


lim (j) L ± 


lim 

t —nitrile 


d4> L 

± dt 


(t>E |an± = 0±> 


lim 

t ^i^max 


dcpE 

r dt 


(38) 

(39) 


In this way by solving the equation of motion in the whole Euclidean domain, the initial 
(or final) data for evolution in the Lorentzian sector can be recovered. It is worth stressing 
that, as already discussed in the previous subsection, the continuity of solutions up to 
at least second time derivatives is required. (The equations are of second order.) The 


additional conditions 


lim 


± 


dt 2 


lim 


d 2 ^ 

t dt 2 


(40) 


therefore have to be satisfied. In general, one can expect that only a subset of all pos¬ 
sible boundary values <f>E |an B obey this condition. Such a possibility might be attractive 
because it may put restrictions on the form of the allowed boundary conditions, based on 
the requirement of mathematical consistency. For the typical equation of motion under 
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consideration, Eq. ©, the requirement (l40]i leads to the following condition of continuity 
of the source term 


lim S[(j) L ±\ = lim S[4>e\ ■ 

t~t—*±t max"^ 


(41) 


In order to illustrate some features of the method discussed above, let us study the 
equation of motion 


d 2 (j) 


m 


d 2 cj) 

dx 2 


0, 


(42) 


with the step-like /3-function 


m 


1 

for 

t t max , 

-1 

for 

Uax U t T t 

1 

for 

t ^ Uax 


(43) 


In this (l+l)-dimensional example the variation of /3 is not continuous across the signature 
change, therefore condition (140|) will have no chance of being fulfilled. Nevertheless, let us 
choose the following boundary conditions: 


(j)± = c± cos(x ) , 


(44) 


where c± G R. are constants. (The boundary condition at <9h2oo will not play any role.) 
With boundary values specified, a solution to Eq. (1-42|) in the Euclidean domain — 
{(t, x) G M 2 ,t max >t> —t max } can be found: 


<Pe = COs(x) 


c+sh(f + t max ) - c_sh(t - f max ) 


sh(2t n 


Using the solution we can find that 


d(f) E , . =Fc T ± c ± ch(2t n 

inn —— = cos x 


t —^ztimax = dt 


sinh(2t ri 


(45) 


(46) 


which will be used to impose Cauchy initial conditions in the Lorentzian domains. Taking 
this into account, solutions to Eq. (T42|) in the domains = {(f, x) G M 2 , ±f > f max } can 
be determined: 


4>l+ 

$L- 


coslx 


cos Ur 


c+ (sh(2f max ) cos (t - t max ) + ch(2t max ) sin(t - t max )) - c_ sin(f - t max ) 


811(24 


(47) 


c_ (sh(2t max ) cos (t + £ max ) - ch(2f max ) sin(f + t max )) + c + sin(t + t max ) 


sh (24 


( 48 ) 


The solutions (j) E and <^l± are presented in Fig. 0 in the form of the density plot for 
the sample values of c± and f max . 
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Figure 5: Density plot of the solutions </>£ and 4>l± to the equation (1421) . The observed 
asymmetry is due to the different boundary conditions at t max and — t max : 4>± = c ± cos(x). 
The values c+ = 1, c_ = 2 and t max = 4 have beed used. 
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The solutions obtained in this way are regular in the whole domain hi = Qg Uf2 + Ufi_. 
However, in agreement with the earlier expectations, the second derivatives are discontin¬ 
uous across the signature change: 


lim 

ax^ 


dt 2 


lim 

t 


d 2 (f) L ± 

dt 2 


2 c± cos(x). 


(49) 


4.2.4 Wave functions 

Mixed-type differential equations for modes are a consequence of properties of effective 
equations in cosmological models of loop quantum gravity. These properties follow from 
a general analysis of consistent and anomaly-free realizations of field equations, so that 
physical predictions made in these models are independent of coordinate or gauge choices. 
The same modifications to space-time structure can be seen from commutator algebras at 
the quantum level, but it is then more difficult to derive dynamical effects. It is nevertheless 
interesting to discuss how signature change or its accompanying instabilities could be seen 
in the evolution of wave functions. 

As long as one has unitary evolution, as can easily be guaranteed in deparameterized 
models which provide the majority of explicit wave-function evolutions in quantum cos¬ 
mology, the wave function cannot be subject to instabilities. (However, we note that at 
a fundamental level it may be difficult to have unitary evolution under all circumstances, 
owing to the problem of time.) Its norm is preserved even throughout high density. Accord¬ 
ingly, the Schrodinger or Klcin-Gordon type equation which the wave function is subject 
to does not change its form to become elliptic; the relevant coefficients in the wave equa¬ 
tion do not change sign. After all, what changes is the signature of space(-time), not the 
signature of (mini)superspace on which the wave function is defined. (See also [SO]. In 
fact, depending on which representation one chooses for the wave function, the equation it 
is subject to may be elliptic or hyperbolic even for standard space-time.) 

Signature change implies instabilities in an initial-value problem for the modes, solving 
effective equations. The fields in effective equations are expectation values of mode oper¬ 
ators in a state obeying the equation for the wave function. At the level of wave-function 
evolution, signature change would therefore be recognized in an exponential change of ex¬ 
pectation values of mode operators. (If wave packets are used, their peak position would 
change exponentially.) For these expectation values (or peak positions), the same sensitiv¬ 
ity to initial values as observed in solutions of effective equations would occur, even if they 
are computed from an evolved wave function. These observable quantities therefore evolve 
in an unstable way, causing the same problem as seen in an initial-value formulation for 
effective equations. Instability problems are just more hidden in schemes that first evolve 
wave functions and then compute expectation values, but they are just as pressing as in 
effective equations. 

In this respect, the situation is reminiscent of quantum chaos, which, compared with 
the classical phenomenon, is more difficult to define and discuss, but not absent, for wave 
functions subject to linear differential equations. In both cases, the sensitivity of evolved 
expectation values to initial choices is relevant. It has been discussed in detail in the 
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context of quantum chaos. Following mm, we can use the same arguments for states 
in the presence of signature change: Unitary quantum evolution ensures that the overlap 
|(0i(f),0 2 (t))| 2 of two different initial states -01 (to) and 02 (to) is conserved in time, and 
therefore does not indicate any sensitivity to choices of initial values. However, quantum 
evolution of classically chaotic systems is very sensitive to perturbations of the Hamiltonian 
operator, in the sense that the overlap |(0 1 (t), 0 2 (t))| 2 changes rapidly if 0i(t) and 0 2 (t) are 
now defined as states evolved from the same initial wave function 0i(to) — 0 2 (to) but with 
different Hamiltonians Hi — H and // 2 = // + eV with some small e and a perturbation 
potential V. Similarly, the same definition of the overlap in Euclidean signature leads to 
rapid (exponential) change because the perturbed evolution by iP 2 contains a factor of 
exp (efrr l Vt). In discussions of quantum chaos, eV is usually thought of as coming from 
unknown interactions with an environment hard to control. The same source of eV exists 
in quantum cosmology because the precise degrees of freedom and interactions at high 
density are not well known. In addition to this, loop quantum gravity is subject to a large 
set of quantization ambiguities, so that the dynamics is not precisely determined and a 
second source of eV results. 

Signature change in effective equations, in all existing models, is implied by the re¬ 
quirement of anomaly freedom, so that predictions are guaranteed to be independent of 
gauge choices. Most models in which one evolves wave functions are based on deparam¬ 
eterization, in which one formulates “evolution” of a wave packet in terms of a so-called 
internal time, which is not a coordinate but one of the dynamical fields. The conditions 
on well-defined evolution after deparameterization are less severe than the general covari¬ 
ance conditions when one does not select a specific time, be it a coordinate or internal. 
If one were to fix the gauge and work in one specihc set of coordinates, one could avoid 
having to use effective equations subject to signature change. Similarly, the possibility 
of formulating stable wave-function dynamics in one chosen internal time does not mean 
that wave functions in general evolve in a stable manner. One would have to show, first, 
that predictions of one’s model do not depend on which field is chosen as internal time 
(a task which has not been performed in deparameterized models of quantum cosmology; 
see [13111]). After this, one could reliably analyze evolution and test whether it is always 
meaningful or has to be stopped when a Euclidean regime is reached. Performing this 
task turns out to be much more complicated than analyzing anomaly freedom of effective 
equations. But interestingly, even in the absence of such an analysis there are hints that 
wave-function evolution becomes unstable at high density, in the sense that expectation 
values of modes change exponentially [ 15] . 

4.2.5 Cosmological implications 

A well-posed treatment of initial and boundary values in loop quantum cosmology implies 
significant departures from the scenarios commonly made in this setting. For instance, 
there has been some interest in a “super-inflationary” phase at high density, around a 
bounce, during which the Hubble parameter grows quickly [H?]. Although the rapid change 
happens for background variables and is therefore not the same as the instabilities of mode 
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equations such as ©, one can easily check in explicit models that the super-inflat ionary 
phase falls within the elliptic range of the partial differential equations for inhomogeneity 
[5]. The rapid growth of background variables cannot be used for cosmological effects 
because their values at the border of the Euclidean phase must be prescribed as part of 
the boundary-value problem. Obtaining a large parameter with an ill-posed initial-value 
problem does not have physical significance. 

A related effect has been seen in a combination of loop-modified background equations 
with inflationary models: At high density, the modified equation for the inflaton has an 
anti-friction term which can easily push the inflaton up to a high value in its potential 
HHIIE]. Given the right potential for suitable inflation, the correct initial conditions can 
therefore be provided. Also this effect is subject to the verdict of being based on ill- 
posed data. An alteration may, however, still be realized if one can show that for a given 
background value 0 at the beginning of the expansion phase there must generically be a 
large 0' according to the well-posed characteristic formulation. Values of 0 far from the 
minimum of a potential could then be achieved afterwards, using well-posed evolution with 
the predicted initial values in the expansion phase. 

With the discussion of initial and boundary values in Sec. 14.2.21 we can make our 
expectations for cosmological scenarios more precise. Some data must be chosen at the 
beginning of the expansion phase, even if the energy density or other coefficients of the 
differential equation never become infinite. This feature is shared with the singular big- 
bang model, in which the main conceptual problem is caused not so much by divergences 
but rather by the presence of unrestricted initial data. (In general relativity, it is possible 
to extend solutions across singularities in a distributional sense, but the extension is not 
unique [49].) If there would be a unique way of deriving initial values at the beginning of the 
expansion phase, the singularity of the traditional big-bang model would not be so much of 
a problem. One could make clear predictions about the initial state and further evolution 
of the universe, as attempted in inflationary scenarios with their (somewhat controversial) 
initial conditions for the inflaton. In practical terms, the singularity problem is therefore 
one of indeterminedness, which is implied but not necessarily equivalent to the occurrence 
of divergences. 

An example for the latter part of the preceding statement is given by signature-change 
cosmology. There are no divergences in these models, and yet an important part of ini¬ 
tial values for the expansion phase remains free. The resulting scenario, which we call a 
finite beginning to emphasize the absence of divergences together with the requirement of 
choosing new initial values, is rather close to the standard big-bang model. In particular, 
for practical purposes the need for new initial values once the density is sufficiently low to 
trigger signature change back to a Lorentzian structure makes the Euclidean region appear 
as a singularity. We therefore speak of a finite beginning instead of a non-singular one, 
keeping in mind that the main problem of a singularity is the indeterminedness of initial 
data. 

There may yet be a crucial difference with the standard big-bang model: What we need 
to prescribe at the finite beginning, according to Fig. [21 is only the field 0, not its normal 
derivative 0' if we start evolution in the collapse phase. If one has some means to know 0 
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in a neighborhood of the border between the Euclidean region and the expansion phase, for 
instance from some hypothetical observations, one can derive q¥ and draw conclusions about 
the field in the collapse phase. There is therefore some connection between collapse and 
expansion, although it is not causal and not fully deterministic. A thorough cosmological 
analysis of scalar and tensor power spectra in signature-change models is required before 
one can tell how much about the collapse phase could be deciphered in this indirect way. 

4.2.6 Global issues 

While effective methods for constrained systems provide reliable local equations, as con¬ 
firmed by the results of this paper, there may seem to be several global problems related 
to the new type of partial differential equations which imply the disappearance of time 
or causality at high density^ One example has been discussed in ra in the context of 
black holes, where the non-singular beginning in cosmology as detailed here is replaced by 
a naked singularity (again in the sense of indeterminism) with a Cauchy horizon. Even if 
local field equations are regular, the lack of a causal structure in some regions may lead to 
unacceptable indeterministic behavior at a global level. 

A different kind of global problem is indicated by one feature of solutions to Tricomi’s 
problem which we have not mentioned so far. Again, locally the equation and its solutions 
are well-defined. But generic solutions turn out to have a pole at one point A D B in 
Fig. Ql where the Euclidean arc ends [7j: Although solutions are smooth in the interior 
of the characteristic region (for smooth boundary data) most of them have a pole at one 
endpoint of the boundary. (This result explains the ubiquity of sonic booms in analogous 
acoustic models.) The solution may remain finite, but not all of its derivatives do. Derived 
quantities such as the energy density in a matter field could therefore diverge. If this 
happens, it is not likely that cosmological perturbation theory gives reliable results, even 
though the local perturbative equations are regular. 

The acoustic analog illustrates a further point: One could think that signature change 
in cosmology (or black-hole physics) is harmless because analogous effects can appear in 
well-known systems such as transonic motion (or the model discussed in Sec. 13. 2\i . The 
difference is that transonic motion leads to signature change only for excitations in the 
fluid, while all other propagating degrees of freedom including the bulk fluid motion and 
space-time are still governed by deterministic equations. (The speed of light is an absolute 
limit for any causal motion, very much unlike the speed of sound in a fluid. If B < 0 in fTTTl . 
all motion is eliminated.) If the bulk of the fluid is forced to move faster than its own speed 
of sound, it overtakes any sound wave generated in it, so that its density prohle provides 
future conditions for the wave. Mathematically, this is represented by Tricomi’s future data. 
However, the bulk fluid itself moves in a deterministic (but not wave-like) way even if it is 
faster than its own speed of sound, and no causality issues appear. The situation is very 
different when signature change happens for space-time physics, in which case no reference 
time remains to define a causal structure and no mode can evolve deterministically. For 

3 “And well I know it is not right// to seek and stay Time in his flight.” [50] 
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this reason, signature change in effective space-time models of loop quantum gravity has 
much more radical consequences than analogous effects in matter systems, as discussed in 
the present paper as well as l-TT I . 

5 Conclusions 

We have described several fundamental properties of a new scenario, based on signature 
change, that has emerged from spherically symmetric and cosmological models of loop 
quantum gravity in recent years. We emphasize that this scenario cannot be seen in the 
minisuperspace models traditionally studied in loop quantum cosmology [52]. In fact, the 
possibility of signature change casts significant doubt on the viability of minisuperspace 
models of loop quantum cosmlogy because in such models one would not see any sign of 
a disappearing causal structure at high density. Minisuperspace models of loop quantum 
cosmology may be used for some estimates of background properties, but they can no longer 
be considered as reliable sources of stand-alone cosmological scenarios. One always has to 
go beyond homogeneity to make sure that there is a well-defined space-time structure and 
to check whether mode equations remain hyperbolic, or to exclude other exotic effects. 

Signature-change cosmology is therefore a new scenario which crucially relies on in¬ 
homogeneous features of loop quantum gravity (even though, needless to say, it has not 
yet been derived from a full quantization of gravity). Its details, regarding for instance 
power spectra, still have to be worked out, but we believe that the more mathematical 
and conceptual properties discussed in this article already show a large number of inter¬ 
esting features. The main consequence in practical terms is the occurrence of instabilities, 
related to the sensitive dependence of solutions on initial data. In some examples, instabil¬ 
ities may be physical effects which imply rapid change but no inconsistencies. In quantum 
gravity, however, instabilities of the kind encoimtered in the presence of signature change 
are fatal: The theory remains subject to a large number of quantization ambiguities in its 
equations, and not much is known about suitable initial states for quantum space (and not 
just quantum matter). With this inherent vagueness, one cannot afford any instabilities 
that would magnify theoretical uncertainties in a short amount of time, even if this time 
may be as small as the Planck time. No predictions would be possible. We are therefore 
sympathetic with the verdict “In light of the fact that even this ’well-behaved’ signature 
change system predicts its own downfall, it may be prudent to reassess the inclusion of 
signature changing metrics in quantum gravity theories.” [53] obtained after a detailed 
analysis of initial-value problems in a model of classical signature change. Our model is 
provided by effective equations of loop quantum gravity in which signature change appears 
to be generic (and is not included by choice). The verdict of [53] can therefore be applied 
to loop quantum gravity at least to the extent that extreme caution is called for when one 
considers evolution through high density. 

In the presence of signature change, instabilities can be avoided only if one switches 
to a 4-dimensional boundary-value problem at high density, giving up causality. Although 
ambiguities remain in the theory, there are several qualitative effects with interesting irnpli- 
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cations. As one of the main observations in this article, the mixed-type partial differential 
equations for modes in this context strike a nice balance between deterministic cyclic mod¬ 
els and singular big-bang models. There are no divergences, and yet initial data in the 
infinite past do not uniquely determine all of space(-time). For every mode, one must spec¬ 
ify one function at the beginning of the expansion phase even if one has already chosen 
initial values for the contraction phase. Still, the normal derivative of the field is not free 
and may carry subtle but interesting information about the pre-big bang. 
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A Space-time 

We present here a somewhat technical discussion of fundamental properties of space-times 
underlying ([I]). Symmetries and gauge transformations are especially important in this 
context, as well as properties of the canonical formulation of gravity. In our following 
exposition, we also provide more details on the derivation and reliability of equations such 
as (P in effective models of loop quantum gravity. 

A.l Gauge transformations 

In spite of its appearance, equation (P) is covariant in a generalized sense, according to an 
effective (and canonically defined) non-Riemannian structure of quantum space-time. The 
usual Lorentz and Poincare symmetries, under which the classical version of (P) with (3=1 
(and constant a) would be invariant, are realized in canonical gravity as a subalgebra of 
the infinite-dimensional algebra (or rather, algebroid [53]) of deformations of 3-dimensional 
spacelike hypersurfaces in space-time. (See for instance [55) [56] .) These hypersurface 
deformations are gauge transformations of any generally covariant theory, including gravity. 

Hypersurface deformations [8] are more suitable than Poincare transformations in sit¬ 
uations in which no background space-time metric is assumed. They provide the proper 
framework for a discussion of generalized space-time structures as the may be implied by 
canonical quantum gravity. Hypersurface deformations have an infinite set of generators, 
spatial ones given by D[N a } with spatial vector fields N a tangential to the hypersurface, 
and normal ones H[N} with functions N on the hypersurface so that the deformation is by 
an amount Nn a along the unit normal n a . The classical space-time geometry [S3 (as well as 
the canonical form of general relativity [5?]) imply that these generators have commutators 
(or classical Poisson brackets) 

{D[N?],D[N b 2 ]} = D[N(D h N* - A^D Q iVf] 

{H[N},D[N a }} = H[N a D a N] 
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(50) 

(51) 


{HimiHim]} 


- JV 2 D 6 JV,)] 


(52) 


with the induced metric h a b and covariant derivative D a on a hypersurface. 

When one tries to quantize the theory canonically, one should turn the gauge generators 
into operators, so that they still have closed brackets given by commutators. Otherwise, 
the classical gauge transformations would be broken and the quantum theory would not 
be consistent; it would have gauge anomalies. The anomaly problem of canonical quan¬ 
tum gravity is important, but also very difficult and unresolved so far. (One reason is 
the presence of the metric in ([52]) . which at the quantum level would be an operator and 
give rise to complicated ordering questions with a quantized D[N a ].) Nevertheless, there 
have been several independent indications in recent years which give some hope that the 
problem can be solved. At the operator level, especially 2 + 1-dimensional models have 
been analyzed in quite some detail, paying attention to the full algebra of gauge genera¬ 
tors. Consistent versions have been found in different ways [5811591 1601 151] . including also 
spherically symmetric models [62] . 

Independently, effective calculations start from the observation that the quantum oper¬ 
ation of commutators together with a closed algebra of operators Cj, [ Cj , Cj) = ffjCx with 
structure constants (or functions/operators) fjj, implies a closed algebra under Poisson 
brackets of effective constraints (Cj). These effective constraints are defined as expectation 
values (Cj) of the constraint operators in an arbitrary state [021 E3J. Effective constraints 
are therefore functions on the quantum state space, which can conveniently be parame¬ 
terized by expectation values and moments with respect to a set of basic operators. In 
addition to these effective constraints, obtained as direct expectation values of constraint 
operators, there is an infinite set of independent ones, derived from the same constraint 
operators as (polCj) for all polynomials in basic operators. All these functions on the 
space of states would be zero on the subspace annihilated by the constraints Cj, imposed 
following Dirac’s prescription. The need for an infinite set of effective constraints for every 
single constraint operator follows from the requirement that a whole tower of infinitely 
many moments must be constrained together with every constrained expectation value. 

A Poisson bracket on the quantum state space of expectation values and moments 
can be defined by {(A), (B)} = {[A,B])/ih, extended by the Leibniz rule to products of 
expectation values (as they appear in quantum fluctuations and higher moments). With 
these definitions, it follows that effective constraints form a closed algebra under Poisson 
brackets if the constraint operators form a closed algebra under commutators. Practically, 
it is easier to evaluate Poisson brackets than commutators, an observation on which the 
idea of canonical effective equations [H3 ES] and constraints [621 E3] is based. Also, there 
are useful approximation schemes, such as a semiclassical one in which one would do 
calculations order by order in the moments, which are easier to implement for effective 
constraints than for constraint operators and allow one to handle the unwieldy space of 
all states more efficiently. For finite orders in the moments, there is a finite number of 
independent effective constraints, and the task of computing their algebra becomes feasible. 

If we have a closed algebra [C/, Cj\ of operators, the corresponding effective constraints 
truncated at some moment order form a closed Poisson-bracket algebra up to this order, 


irrespective of the states used. This formulation of effective theories is therefore much more 
general than the usual idea of effective actions [67] . (The latter are often combined with 
further approximations such as derivative expansions, or with restrictions of states such as 
near-vacuum states for the low-energy effective action.) In our case, we need not make any 
assumption on the class of states in order to obtain a closed algebra of effective constraints. 
One may just have to include higher orders in the moments for a better approximation to 
solutions corresponding to strong quantum states; but even at lower orders, the algebra 
must be closed (to within the same order). Effective constraints therefore provide a good 
test of possible anomaly-free quantum constraint algebras. In particular, if one can show 
that no closed effective constraint algebra exists to within some order in moments, using 
a large-enough parameterization of quantum corrections, there cannot be a closed algebra 
of constraint operators. And if closure of effective constraints can be achieved only if the 
classical constraint algebra is modified, the full quantum constraint algebra must be subject 
to quantum corrections which change the form of gauge transformations. Signature change 
is just one of the consequences found in this way: Instead of (152]) . we then have 

{HlN^HlNz}} = —D\fih ab (NiD b N 2 - N 2 D b N 1 )} (53) 

with a phase-space function /3 which may turn negative (while the other brackets involving 
D[N a ] remain unchanged). Canonical held equations consistent with this modified bracket 
are of the form (CQ), as shown by [5] using the methods of [8; ;68] • (The same modification 
of derivative terms is realized to higher orders in derivatives [69].) 

So far, calculations in models of loop quantum gravity have been performed only to 
zeroth order in the moments. But they still show interesting effects because, in addition to 
moment terms, the theory implies further quantum corrections. At high density, holonomy 
modifications are relevant, which are implied by the basic assumption of loop quantum 
gravity m mm that the gravitational connection can be represented as an operator 
only when it is integrated and exponentiated to a holonomy. A second effect, inverse-triad 
corrections [SI USE], is more indirect but also related to the basic assumption. It implies 
corrections relevant at lower curvature PISIIS]. Effective constraint algebras have been 
computed in both cases and found to allow for consistent versions: [SI [SUS! as examples 
for inverse-triad corrections, [3] for holonomy corrections, and ID 1212] for combinations 
of both. Whenever holonomy modifications are present, all generic consistent effective 
constraints found so far imply mode equations of the form ([Tj). 

The derivation shows that covariance is not broken because all classical gauge gener¬ 
ators have a valid analog as a constraint operator Cj or as an effective constraint. For 
every classical gauge transformation, there is a corresponding quantum or effective gauge 
transformation. No transformations are violated, including the Lorentz and Poincare ones 
that one obtains as a special case of hypersurface deformations [55] . Nevertheless, it is 
possible for (JT]) to differ from standard covariant wave equations because quantum gauge 
generators, while they must not be broken for an anomaly-free theory, may be subject to 
quantum corrections. These corrections can be seen in the structure operators or func¬ 
tions ffj of the quantum or effective constraint algebras. For holonomy and inverse-triad 
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corrections, the classical structure functions are, as in (153]) . multiplied with a phase-space 
function /?, which determines a new consistent form of quantum space-time covariance. 
The same function appears in mode equations (JT]) derived from the effective constraints. 

A.2 Slicing independence 

Even though the speed in (jT]) provided by quantum geometry depends on the spatial metric 
and extrinsic curvature — quantities that are not covariant in classical space-time — the 
effect is frame independent in a subtle way. The models in which such modified speeds have 
been derived are anomaly free, so that the classical set of gauge transformations (given by 
coordinate changes) is not violated. However, not just the dynamics but also the structure 
of space-time receives quantum corrections. There is a new set of gauge transformations 
under which the quantum-corrected held equations are invariant, and which has the full 
set of standard coordinate changes as its classical limit. The effective theories, including 
modified speeds they predict, are covariant under these deformed gauge transformations. 
Accordingly, the space-time structure is no longer classical or Riemannian, but it remains 
well-defined in canonical terms. 

Put differently, one could worry that models in which waves propagate with speeds 
depending on the spatial metric and extrinsic curvature of constant-time hypersurfaces 
violate the slicing independence of the classical theory. Classically, for any given space- 
time, one can, depending on one’s set of coordinates, choose equal-time slices with large 
or small extrinsic curvature, even if the covariant curvature of space-time vanishes. A 
speed depending on extrinsic curvature would then suggest that predictions depend on the 
slicing or the choice of an initial-value surface, which would be unacceptable. This worry 
is unjustified because the two initial-value surfaces, one with small and one with large 
extrinsic curvature, would give rise to different physical solutions of the effective theory 
even though they would present the same classical solution. Predictions, for instance of 
the speed, then differ simply because one would consider physically different solutions. 
Transformations between different slicings are, just like coordinate changes, gauge trans¬ 
formations which are modified in the effective theory. Slicings of one and the same classical 
space-time, which are related by classical gauge transformations, are not gauge related in 
the quantum-corrected setting. 
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